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In quantum communication and distributed quantum computing, one-dimensional waveguides provide direc- 
tional transfer of quantum information. A single-mode waveguide has a density-of-states singularity at the lower 
cut-off frequency, which resembles sharp resonances of a cavity but with non-Markovian dynamics. Thus we put 
forward schemes of coupling atomic transitions and a waveguide continuum edge. We first present a scheme of 
spin-photon quantum interface operating in the non-Markovian regime for a A-type three-level system coupled 
to a waveguide. Then we show that strong coupling between atomic transitions and a waveguide continuum 
edge can lead to vacuum Rabi oscillations and bound polariton states. 
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I. INTRODUCTION 

Quantum networks are essential in quantum communica- 
tion and distributed quantum computing fUQ. A quan- 
tum network consists of local nodes and connecting quantum 
channels. The stationary qubits in local nodes can be provided 
by trapped ions ISJCI. collective atomic excitations 151 121, 
quantum dots ID [10] [TT], or solid-state impurities lfT2HT6l . 
Single photons produced by manipulation of atomic transi- 
tions in such systems can form the flying qubits. The quan- 
tum information carried by the flying qubits can be conducted 
between local nodes through waveguides. To improve the ef- 
ficiency of conversion between stationary and flying qubits, 
cavities have been adopted to enhance the photon emission |2 - 

aiiiiEi. 

In this paper, we present schemes of strong coupling be- 
tween atomic transitions and waveguide continua without an 
intermediate cavity. The basic idea is: at the lower cut-off 
frequency of a single-mode waveguide, there is a singularity 
in the density-of-states (DOS) of the one-dimensional con- 
tinuum. Such a DOS singularity resembles the sharp reso- 
nances associated with discrete states of cavities. The line- 
shape of the singularity, in contrast to those of cavity res- 
onances, is highly non-Lorentzian, and therefore the photon 
emission at the continuum edge is a non-Markovian process. 
Such non-Markovian nature induces interesting memory ef- 
fects on, e.g., spin-photon quantum interfacing. We derive 
an exact solution for the quantum interfacing scheme in the 
non-Markovian regime, extending the control scheme of spin- 
photon interfaces that works in the Markovian regime ifTTI . 
Also, the coupling between atomic transitions and waveguide 
continuum edges can lead to vacuum Rabi oscillations and 
bound polariton states, similar to atom-cavity quantum elec- 
trodynamics in the strong coupling regime. These results are 
relevant to recent investigation on coherent coupling between 
atomic transitions and one-dimensional continuum of photons 
orplasmonics lfT8l420l . 



'Electronic address: |rbUu@phy.cuhk.edu.hk| 



There have been a number of proposals |l2l [3] |8] |9] [TTII and 
experimental implementations f2P-54l of quantum interface 
between stationary qubits and flying photon qubits. A mile- 
stone is the proposal of the state transfer between two atoms 
in cavities |2|. The scheme is based on the idea of using two 
mutually time-symmetric laser pulses to control the cavity- 
assisted Raman processes in the two atoms at the sending and 
receiving nodes. Such a scheme, however, requires the two 
nodes be identical and is sensitive to randomness in atom- 
cavity couplings. In order to overcome such obstacles, adi- 
abatic schemes similar to the stimulated Raman adiabatic pas- 
sage ||25l have been proposed to implement the cavity-assisted 
Raman processes |[3][8l|9|. The adiabatic scheme has the ad- 
vantage of robustness against parameter uncertainty but the 
disadvantage of slow operation rates. An exact solution of 
the control pulse for spin-photon interfacing was later discov- 
ered IfTTI . in which the photon wavepacket can be arbitrarily 
tailored by shaping the laser pulse that controls the Raman 
process in a A-type atomic system coupled to the cavity. By 
combining the sending and receiving functions of two spin- 
photon interfaces, deterministic state transfer and entangle- 
ment between two remote stationary qubits are possible, with- 
out requiring two identical nodes. Such a control scheme is 
particularly suitable for solid-state quantum networks, where 
local nodes can hardly have identical parameters. In the above 
spin-photon interface schemes, the cavity modes are coupled 
to a wide band of photon continuum in free space or waveg- 
uides. Thus the photon sending or receiving can be well de- 
scribed by the input-output theory |26| with Markovian ap- 
proximation (i.e., Wigner-Weisskopf approximation) |27|. 

The quantum interface scheme proposed in this paper uses 
direct coupling between atomic transitions and the DOS sin- 
gularity of a waveguide continuum (without a cavity as the 
intermediate). The photon emission/absorption becomes a 
non-Markovian process due to the non-Lorentzian lineshape 
of the photon continuum DOS. We extend the exact solution 
in Ref. IfTTI to the general case of non-Markovian coupling. 
The exact solution reduces to the input-output theory in the 
limit of Markovian approximation. In order to investigate the 
validity of the Markovian approximation, we numerically sim- 
ulate the control scheme for various operations for different 



coupling parameters (such as detuning frequencies, coupling 
constants, and pulse durations). The quantum interfaces in- 
cluding waveguides are widely used in solid-state systems, 
such as cooper pair boxes coupled to a transmission line Ii28l 
and single emitters coupled to the surface-plasmon modes of a 
metal nanowire ifTSl . Quantum dots can also couple to etched 
waveguides or ID photonic crystals |'l9"29l to form quantum 
interfaces. Our scheme works for those quantum interfaces in 
the strong coupling regime, where the Markovian approxima- 
tion may not be valid. 

The DOS singularity of a waveguide continuum resembles 
a sharp cavity resonance. Therefore we envisage strong cou- 
pling between atomic transitions and the continuum edge. We 
study the photon emission dynamics of a two-level system, 
and find vacuum Rabi oscillations of the non-Markovian emis- 
sion and a bound polariton state which results in photon local- 
ization and incomplete emission. These results are related to 
the photon localization phenomena previously discovered in 
photonic band-gap materials |28, 30, 31 1. 

This paper is organized as follows: In Sec. [ll] we explain 
the basic idea of quantum interface in the non-Markovian cou- 
pling regime, present an exact solution of the control, and nu- 
merically investigate the validity of the Markovian approxi- 
mation and the effects of photon leakage. In Sec. [Ill] we study 
the spontaneous emission of an exited two-level system cou- 
pled to a one-dimensional waveguide. Sec. IV 
and concludes the paper 



summarizes 



II. SPIN-PHOTON INTERFACE BY ATOM- WAVEGUIDE 
QUANTUM ELECTRODYNAMICS 

A. Description of the control sclieme 
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(c) 




(b) 



FIG. 1: Illustration of atom-waveguide quantum electrodynamics, 
(a) The local nodes composed of three-level systems. |1) and 
|2) (|2'» are two-near degenerate ground states forming the qubit and 
|3) (|3'» is an excited state. The state |2) (|2')) is coupled to the inter- 
mediate state |3> (|3'>) by the waveguide modes {\k)] with strength gj. 
{g'f,) and the state |1> (|1'» is coupled to the state |3) (|3'» by a classi- 
cal pulse with Rabi frequency Sl(t) [Q.'(t)]. These two nodes are con- 
nected by the waveguide. The output of the sending node is directed 
to the receiving node as its input, (b) The two-level system embed- 
ded in a one-dimensional waveguide. The transition couples to the 
waveguide modes with strength g^. (c) Schematic density of states in 
a one-dimensional waveguide near the lower cut-off frequency (iJo). 



Figure [T] (a) shows the scheme of quantum interface and 
state transfer. The basic idea is to control the emission of a 
single photon in the sending node into the directional waveg- 
uide and the absorption of the photon at the receiving node so 
that quantum information is transferred between the two sta- 
tionary qubits The stationary qubit in a node is encoded in 
the two near degenerate ground states of a A-type three-level 
system (|1) and |2) in the sending node, and 11') and |2') in the 
receiving node). The two three-level systems are connected 
by a quantum channel formed by a semi-infinite waveguide. 
The coupling to the photon continuum in the waveguide is 
mediated by the excited states |3) and |3') in the sending and 
receiving nodes, respectively. Note that the two nodes are not 
required to be identical. The driving laser with Rabi frequency 
Q. is resonant to the transition |1) <-> |3). If the sending node 
is initially in the state |1), the laser pulse is designed such 
that the system is excited to the intermediate state |3) and then 
relaxes to the state |2) with a photon wavepacket emitted to 
the waveguide. The photon wavepacket can be shaped on de- 
mand by designing the control laser pulse Q(f). When the 
photon arrives at the receiving node, a control laser pulse can 
be designed to absorb the photon without reflection (quantum 
impedance match). The design of the control pulse at the re- 
ceiving node for a certain arriving photon wavepacket can be 



obtained by time-reversal of a control laser pulse for sending 
the time-reversed photon wavepacket IfTTI . The DOS singu- 
larity of the photon continuum in the ID waveguide provides 
strong coupling to facilitate the quantum interface (so that the 
leakage of quantum information into free space is negligible), 
which otherwise requires a cavity to enhance the coupling. 
Near the continuum edge of the ID waveguide, the photon 
emission and absorption are in general non-Markovian. A key 
issue to be addressed in this paper is how to design the con- 
trol laser pulse Q(f) to realize the emission or absorption of 
an arbitrary photon wavepacket in the non-Markovian regime, 
which is the inverse problem of quantum evolution of the cou- 
pled atom-waveguide system under time-dependent control. 



B. Exact solution for control design 

In this subsection, we present the exact solution of the con- 
trol laser pulse for an arbitrary photon wavepacket. Here we 
neglect for the moment photon leakage into free space through 
the excited states |3) and |3'). The control errors induced by 



the free-space emission will be investigated later in Sec. II D 



The quantum interface is shown in Fig. [TJa). The state |1) 



3 



is coupled to |3) with x -polarization light. And the state |2) 
is coupled to |3) with y -polarization waveguide modes. The 
Hamiltonian describing the interaction between the three-level 
systems and the waveguide continuum, with rotating-wave ap- 
proximation and in the rotating reference frame, is 



H = 



k 

J]g^|3')<2'|fl,e'('=32-^)' + Q'(f)|3')<l'| 



+ h.c. 



+ h.c., (1) 



where at is the annihilation operator of the waveguide mode 
\k), the transition |2) «-> |3) is coupled to the waveguide mode 
\k) with strength gk, the transition |1) «-> |3) is driven by 
a laser pulse with time-dependent Rabi frequency Q.(f) with 
central frequency equal to the transition frequency between 
|3) and |1), £32 is the energy splitting between |3) and |2). 

LOk — -^Wq + {c/n)^k^ is the photon frequency with loq denot- 
ing the lower cut-off frequency, and n denoting the refractive 
index of the waveguide. Near the cut-off frequency, the pho- 
ton dispersion is close to that of a massive particle and the 

DOS = ^(?7;i) ^= has a singularity [Fig.jljc)]. Symbols 

with prime marks have similar meaning but for the receiving 
node. 

The coupled system has an invariant Hilbert sub- 
space expanded by the basis states {|l,2')|vac), |3,2')|vac), 
|2, 2')\k), |2, 3')|vac), |2, r>|vac». Here the front vectors de- 
note the states of the three-level systems, |vac) denotes the 
vacuum state of the waveguide, and \k) stands for the one- 
photon Fock state of the waveguide mode of wavevector k. 
The state-transfer works for arbitrary initial state of the send- 
ing node, while the receiving node should be initialized to the 
state |2'). The state of the system can be written as 



mt)) = a\2,2')+mt)). 



(2) 



with 



|^(f))=Ci(f)|l,2'> + C3(f)|3,2'>-H C2(k)\2,2'mdk 

Jo 

+ D3(0|2,3'>+Di(f)|2,l'>. (3) 

The coefficients satisfy the Schrodinger equation 

d,Ci(t) = -/Q*(0C3(0, (4a) 

d,C3(t) = -/Q(f)Ci(f) - / gkC2(k)e'''"-"^'''dk, (4b) 
Jo 

d,C2(k, t) = -iglC3(t)e''^"-''''-' - ig[*D,(t)e'"''-''y, (4c) 

Jr»+oo 
g'^Czikye'^y-'^^'dk, (4d) 




5,Di(f) = -/Q'*(0£»3(f). 



(4e) 



If the waveguide is long, we can separate the state trans- 
fer into two independent steps, namely, sending and receiving 
of the photon wavepacket. The emitted photon at the remote 
future for the sending node can be treated as the incoming 



photon at the remote past for the receiving node. The output 
wavepacket of the first node constitutes the input one for the 
second node with a time delay. We regard the time delay to 
as ^, where L is the length of the waveguide and Vg is the 
group velocity of the photon wavepacket in the waveguide. 
The waveguide is long enough so that the second atom is de- 
coupled from the waveguide during the photon emission in 
the sending node. When the wavepacket propagates down the 
waveguide and enters the second node, the sending node does 
not interact with the wavepacket any more. Therefore, the dy- 
namic equations are decoupled into two parts for the sending 
and receiving process respectively, 

5,Ci(f) = -!0*(f)C3(f), (5a) 

/-*+oo 

d,C3(t)^-imCi(t)-i g,C2(k)e"''^-'-"'''dk, (5b) 
Jo 



d,C2(k,t)^-iglC,(t)e'"^'-'''''-'. 
for f a; , and 



(5c) 



d,Di(t)^-in'*(t)Di(t), (6a) 

r*+oo 

d,D3(t)^-i£l'(t)Di(t)-i g',C2(k)e'''^^-'-'''''dk, (6b) 
Jo 

d,C2(k, f) = - ig';D,(t)e"^^'-''y, (6c) 
for f X t(). 

For the sending node, we look for a solution of the driving 
pule Q.{t) for an arbitrary output F(k), which determines the 
photon wavepacket shape. The boundary conditions are taken 
as C2{k, -00) = and C2ik, +00) - F(k). This means that the 
sending node functions under the condition of no incoming 
photon and generates an outgoing photon wavepacket of the 



desired shape. From Eq. (5c 1, we have 



F{k) 



•J —cx 



C3(f)^^e'"''-''^^-'df. 



(7) 



By inverse Fourier transform, the coefficient C^it) is obtained 
as 



C3(0 



Jo gv 



F(k)e-'""^'''^-' — , 
Si 27r 



(8) 



from which C2(k, t) is derived from Eq. (5c 1. To get the coef- 
ficient C\{t), we write it as |Ci(f)|e'*'*'', where the phase (f)\ is 
a real number. The amplitude of Ci(f) is given by the normal- 
ization condition |Ci(f)l^ = 1 - |C3(f)l^ - \C2(k, f)\^dk and 
the phase 



\J —CO 



, |C3(f')l'5,03 - /o" \C2{k, n\^dr<P2dk 



>\(t)= df r . (9) 

Finally, from Eq. (j5]), we derive the driving field Q(?) in terms 
of all coefficients. 



fi(f) = [ic-^\t)d,c,(t)\ . 



(10) 
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The control of the receiving node is similar to that of the 
sending node but with different boundary conditions. The 
photon emitted by the sending node is fully absorbed by 
the receiving node. This means the boundary conditions 
C2(k, to + t') = F{k) for f -> -oo and C2(k, to + f') = for 
t' +00 (here for convenience we have written t = to + f ). 
The coupling strength between the atom transitions and the 
quantum channel depends on the position of the node through 
the relation Ell ^ = g*e'*^. The coefficient Di{t) and Ctit) 
is expressed as 

Diit' + to) = -/ —Fik)e-'^"'-''^2->«'^'«) dkaJk, (11a) 
Jo gt 2;r 



C2(kj'+to)^F(k)- 



f dt" f 

J-co Jo 



dk' 



In 



For the receiving process, the control pulse is obtained as 

Q'(f' + fo) = [iD-,\t' + tQ)d,Di(t' + to)]* . (12) 

By combing the sending and receiving processes, the quantum 
state transfer can be implemented and remote entanglement 
can also be generated ||TT|. Thus the exact solutions give a 
general control scheme for quantum interface. 

In the Markovian approximation, the equations of motion 
are decoupled into the following equations 11261 



(13a) 
(13b) 



d,Ci(t) = -in*(t)C3(t), 

drCi(t) = -iD.(t)Ci(t) - |C3(f), 



where the decay rate y = 27r|gi.|^(^) is a constant, in which 
the DOS (dk/dcok) is taken as the value at the central fre- 
quency of the output wavepacket. The approximation relies 
on the assumption that the DOS of the waveguide is flat within 
the outgoing wavepacket spectral width. In the Markovian ap- 
proximation, the control pulse Q.{t) is still given by E q. ([T0| , 
but the coefficients (Ci and C3) are determined by Eq. (|13|i. 



C. Numerical investigation 

In this subsection, we numerically study the control design 
of quantum interface. In particular, by numerical comparison 
of the control designs with and without the Markovian ap- 
proximation, we investigate the validity range of the approx- 
imation. The results show that the Markovian approximation 
is well justified only when the spectral width of the photon 
wavepacket is much less than the separation of the wavepacket 
center frequency from the waveguide cut-off frequency. For 
the moment, the photon leakage into the free space is assumed 



to be negligible, but will be considered later in Sec. 11 D 



In the simulation, we assume that the photon wavepacket 
have the form 



Fik) = C X sech 



(14) 



where C is the normalization factor to ensure that 
/ \Fik)\^dk = 1, a>i is the central frequency, and ctq is the 
spectral width of the photon wavepacket. To be specific, we 
take the model system as a doped quantum dot coupled to a 
rectangle waveguide [ 1 1 1 . The two ground states |1) and |2) 
are represented by the two electron spin states. The excited 
state 1 3) is a trion state formed by two electrons in a singlet 
state plus a heavy hole. The dipole moment which decides the 
coupling strength is assumed to be 75 Debye and the cut-off 
frequency of the waveguide is set to be 1.5 eV. We consider 
the sending node with the initial state in |1). 

Figure |2] compares the control schemes obtained with and 
without the Markovian approximation. The central frequency 
of the outgoing wavepacket is chosen to be 1 meV above the 
waveguide cut-off frequency, and £32 = 1.501 eV. The spec- 
tral width of the wavepacket is 0.08 meV. Fig.|2|a) compares 
the control pulses obtained with and without the Markovian 
approximation. At the beginning stage of the control pro- 
cess, the two control pulses are almost the same. In the later 
stage, the exactly obtained control pulse lasts longer than the 
Markovian approximation to complete the full Raman transi- 
tion. This can be understood by memory effects in the photon 
emission due to the DOS singularity structure in the waveg- 
uide continuum. Such memory effects behave as if the emitted 
photon can be re-absorbed due to the reflection by the waveg- 
uide walls. Thus even after the photon emission, a strong 
laser pulse is still required to avoid the Raman process be re- 
versed. Fig. |2](b) shows the state transfer calculated without 
the Markovian approximation. After about 20 ps, the Raman 
process for the photon emission is almost completed and the 
population of the state |2) approaches one. Fig. |2] (c) illus- 
trates the output wavepacket calculated with the exact equa- 
tions of motion [Eq. (|5]l] but the control pulse designed us- 
ing the Markovian approximation. A small error in photon 
wavepacket generation is induced as demonstrated by the ap- 
pearance of a non-zero imaginary part of the wavepacket func- 
tion. 

Figure. [3]shows the control under the same conditions as in 
Fig. |2]but with the spectral width of the photon wavepackt 
to be much smaller (ctq = 0.008 meV). In this case, the 
photon emission is much slower than in Fig. [2] since the 
wavepacket duration is inversely proportional to the spectral 
width. Indeed, Fig. |3] (b) shows that the occupation of the 
exited state during the Raman process is very small, which 
indicates that the process is nearly adiabatic. In such a slow 
process, the Markovian (memory) effect should be less im- 
portant and the photon re-absorption is negligible. Thus after 
the photon emission, the control is insensitive to the driving 
pulse, which just needs to be turned off. As shown in Fig. [3] 
(a), the Markovian approximation produces almost the same 
control pulse as the exact solution. The validity of the Marko- 
vian approximation in this case can also be understood by 
the fact that the DOS of the waveguide continuum is nearly 
flat within the spectral width of the photon wavepacket even 
though the central frequency is quite close to the cut-off fre- 
quency ((To « wi - <jJo)- This is seen in Fig. [3] (c) which 
plots the DOS against the spectrum of the photon wavepack- 
ets obtained with and without the Markovian approximation. 
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FIG. 2: (color online) Control of a sending node to emit a photon 
wavepacket. The photon wavepacket is 1 meV above the cut-off 
frequency of the waveguide and the spectral width is 0.08 meV. (a) 
The Rabi frequency of the driving pulse designed with and without 
Markovian approximation, plot in black solid and red dashed lines, 
respectively, (b) Populations in the ground state 1 1> (blue dotted line), 
the ground state |2) (red dashed line), and the exited state |3) (black 
solid line), calculated without the Markovian approximation, (c) The 
outgoing photon wavepacket. The red dashed line is the ideal photon 
wavepacket (with zero imaginary part), and the black solid lines with 
symbols are the outgoing wavepacket obtained with the Markovian 
approximation, which has a non-zero imaginary part (black line with 
triangle symbols). The blue solid line is the DOS in the waveguide. 



Fig. [3] (c) also confirms that the photon wavepacket derived 
with the Markovian approximation is almost the same as the 
exact result under the slow emission condition. 



D. Effect of photon leakage 

Now we consider the control errors due to photon leakage 
into the free space. The leakage reduces the waveguide cou- 
pling efficiency. This loss is taken into account by adding a 
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FIG. 3: (color online) The same as Fig. |2] but the desired photon 
wavepacket has much smaller spectral width (0.008 meV). Here the 
imaginary part of the outgoing wavepacket obtained with the Marko- 
vian approximation is equal to zero. 



decay parameter y/ in the dynamic equations, 



d,Ci{t) = 



- /Q*(f)C3(f), (15a) 
- iD.(t)Ci(t) - ij gkC2(k)e'^'''-''''^'dk 



yr 



C3(f), 



d,C2{k,t) ^ - iglC2{t)e 



ltiJl(t-lE32t 



(15b) 
(15c) 



We define the fidelity of the sending operation as 
\{F(ky'^'"''\F(k))\, where F(ky'^''^^ is the ideal sech pulse, and 
F(k) = C2{k, t — > -Hoo) is calculated by Eq. ( 15 i. The control 
pulse in Eq. ( [T5| is obtained by Eq. ( 10 1 without the Marko- 
vian approximation. 

Table. |l] gives the fidelity for sending a photon wavepacket 
with various spectral widths and leakage rates. The param- 
eters are the same as in Fig. |2] and Fig. [3] unless otherwise 
specified. When the leakage is up to 6 % of the emission rate 
to the waveguide, coiTesponding to numerical calculation in 
Ref. fJSI, the operation fidelity is still above 96%. When the 
quantum dot is put in the free space, the free-space emission 
rate is 1 % of the emission rate to the waveguide, and the fi- 
delity is above 99%. The fidelity is determined by the ratio of 



6 



the decay rate and the leakage rate, insensitive to the spectral 
width. 



TABLE I: Fidelity of the sending operation with finite photon leak- 
age rate (7'). The parameters are the same as in Figs. 2 and 3. For 
such parameters, the emission rate into the waveguide (estimated 
with Markovian approximation) is 7=0.27 meV. 



spectral width (meV) 


0.08 0.08 0.008 0.008 


leakage 7' 


1 % 7 6 % 7 1 % 7 6 % 7 


Fidelity 


0.9916 0.9667 0.9900 0.9606 



III. STRONG COUPLING EFFECT 

When coupled to a continuum, a discrete state becomes a 
resonance with finite lifetime. If the DOS of the continuum is 
flat, the discrete state decays exponentially as a typical Marko- 
vian process. If the DOS of the continuum varies abruptly, 
such as the ID photonic continuum near the cut-ofF frequency, 
the resonance of the discrete state does not have a Lorentzian 
shape. Therefore the non-Markovian effect is important in the 
dynamics of such systems. If the DOS of the continuum ap- 
proaches a 6 function, such as in the case of a discrete cavity 
mode, the coupling leads to the level splitting and the vacuum 
Rabi oscillation. A single-mode waveguide has a singularity 
in the DOS at the lower cut-ofF frequency, which may replace 
a discrete state of a cavity. Thus we envisage a scheme of 
strong coupling between atomic transitions and a waveguide 
continuum edge. In this regime the decay of an exited atom is 
no longer exponential as in the Markovian approximation. 



A. A two-level atom-waveguide model 

The model system is schematically shown in Fig. [T](b). A 
two-level atom is embedded in the waveguide. |1) and |0) are 
the exited and ground states of the atom, respectively. The 
transition |1) <-> |0) is coupled to the waveguide mode \k} with 
strength g. Once the atom is in the exited state, it will sponta- 
neously decay by emitting a photon into the waveguide. With 
the rotating wave approximation, the Hamiltonian of the sys- 
tem is 

H = ^ fiUka^ak H — tiu\[)Cr, + ti ^ gicr+ak + cr^al), 
k ^ k 

where coiq is the transition frequency of the atom, cr, - 1 1 )( 1 1 - 
|0)<0|, 0-+ = |1)(0|, and cr_ = |0)<1|. We assume that the 
coupling strength g is a constant in the following calculation 
for simplicity. Here we do not include the free-space emission 
of the atom and the loss of the waveguide. 



B. Theoretical formalism 

We use the Green's function to study the dynamics of the 
system. The advanced propagator G+ and retarded propagator 



G are expressed as ll34l 



1 



H ± /0+ 



(16) 



We assume that initially the atom is in the excited state 
and the waveguide mode is in vacuum state, i.e., |^o) = 
|l)|vac). The initial state is the eigenstate of the noninteract- 
ing Hamiltonian Hq - Yjk^'^kO^cik + jtnoiocr-. We define 
the interaction part as V = H - Hq. The matrix element 
G{cl>) - {tpQ\G\(fo) is calculated by summing all perturbative 
expansion terms of G{co) in powers of V, 



G(aj) 



1 



a> - a)\Q - R{aj) 



(17a) 



where the self energy R(cd + iO^ ) = A(cli) + i with the imag- 
inary part r{(L)) being the transition rate from the discrete state 
to the continuum, and the real part A(aj) the energy shift of the 
discrete state due to the coupling. The explicit expressions of 
r(w) and A(w) are 134) 



r(w) = Ing^DOSico), 



(18a) 



A(^^) 



-"-^ f r(w') — ^- — dw', (18b) 

ZTI J (Jj - Ci>' 



where ^ means the principal part of the integration. 

The Green's function gives the evolution operator in the fre- 
quency domain. 



t/(w) = -^[G_(w)-G+(w)] 
2m 



1 



n[aj- wio - A(w)]2 + [r(co)/2f 



(19) 



And the Fourier transform gives the the time-dependence of 
the probability amplitude in the excited state. 



Ui(t) 



(20) 



In the weak coupling regime, the denominator in Eq.( 19 1 is 
large except in the neighborhood of wiq. We replace A{ll>) and 
r(w) by A(wi()) = A and r(wio) = F. Therefore, the probabil- 
ity amplitude in the excited state is given by the Weisskopf- 
Wigner approximation. 



-;(a>io+A> 



e 



(21) 



C. spontaneous emission in the waveguide 



Figure |4] compares the spontaneous emission obtained by 
Markovian approximation with the exact results, for various 
atom dipole moments. In the strong coupling regime, the vac- 
uum Rabi oscillation is observed [Fig. |4] (b)], which we as- 
cribe to the splitting between a discrete bound state and a DOS 
singularity of the photon continuum. At the long-time limit, 
the oscillation is damped but some residual population at the 
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And the photon wavepacket /(x, t) is calculated by the Fourier 
transform 



t(ps) 



2 3 
t(ps) 



FIG. 4: Populations at the excited state of a two-level atom with 
a dipole moment p of (a) 75 and (b) 300 Debye. The transition 
frequency is 1 meV above the cut-off frequency (1.5 eV), and the 
Markovian spontaneous emission rates are 0.27 meV in (a) and 
4.37 meV in (b). The red solid lines are the exact results and the black 
dot lines are the Markovian approximation. The blue dashed line in 
(b) is the exact result with the photon leakage rate of 0.033 meV. 





1200 




1000 




800 


C/3 




Q. 




CD 


600 


E 




i- 






400 




200 



2 4 6 8 

Position (mm) 

FIG. 5: (color online) The photon wavepacket amplitude in the 
waveguide at various times after the trigger of the emission at f = 0. 
The atom is placed at x = 0. 



exited state of the atom is still left. Such incomplete emis- 
sion evidences the formation of a bound polariton state. The 
residual atomic population increases as the coupling strength 
increases, on the account of the enhancement of the atomic 
component of the polariton state. 

The formation of a bound polariton state leads to the photon 
localization. To show this effect, we calculate the real-time 
evolution of the photon wavepacket during the emission and 
propagation (Fig. |5]l. We expand the state \(p{t)} of the two- 
level system into the following superposition. 



|^(f)) = Ci(f)|l>|vac>+ d^Co(^,f)|0)|^> 



f(x,t) 



dkCo(k,To)e 



-iutit-To) ikx 



(23) 



(22) 



where Tq is the final time of the emission taken as 10 ps. As 
shown in Fig. |5] the photon wavepacket propagates slowly af- 
ter emission (at a speed only 1 % that of light in vacuum), and 
some portion of the wavepacket is localized at the atom posi- 
tion, which is related to the photon component in the bound 
polariton state. 

When the realistic leakage emission into free space is con- 
sidered, the bound polariton state has a finite lifetime and the 
residual atomic population will eventually decay to zero after 
the photon emission into the waveguide. The numerical calcu- 
lation shows that when the ratio of free-space emission to the 
waveguide emission is less than 1%, the free-space leakage 
makes no significant changes to the strong coupling effects as 
shown in Fig|4](b). 



IV. CONCLUSION 

In summary, we have studied a general control scheme of a 
spin-photon quantum interface consisting of a three-level sys- 
tem coupled to a ID photonic continuum and extended the ex- 
act solution of the control scheme 1 11 1 to the non-Markovian 
regime. We found that the Markovian approximation only ap- 
plies to relatively slow operation. For rapid operation near the 
cut-off frequency of the waveguide, the non-Markovian dy- 
namics is essential. We have also studied the effect of imper- 
fections such as photon leakage and found high fidelity for the 
operation in presence of realistic free-space emission. With 
realistic leakage to the free space of 6%, it is possible to main- 
tain the fidelity above 96% for an operation time ranging from 
100 ps to 800 ps. 

In addition, for a coupled atom-waveguide system, we have 
studied the photon emission of two-level atoms and in partic- 
ular the effect of the DOS singularity at the photon continuum 
edge. We show the vacuum Rabi oscillation and the formation 
of a bound polariton state, which reflect the strong coupling 
between a discrete atomic transition and a photonic contin- 
uum edge. We found the photon localization effect due to the 
formation of a trapped polariton. This effect realizes strong 
atom-photon coupling without a cavity, with the DOS singu- 
larity in a waveguide mimicking a discrete state in a cavity. 
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